arXiv:1503.06324vl [quant-ph] 21 Mar 2015 


Convergence and adiabatic elimination for a 
driven dissipative quantum harmonic oscillator* 

R. Azouit ^ A. Sarlette* R Rouchon^ 

March 24, 2015 


Abstract 

We prove that a harmonic oscillator driven by Lindblad dynamics where 
the typical drive and loss channels are two-photon processes instead of single¬ 
photon ones, converges to a protected subspace spanned by two coherent 
states of opposite amplitude. We then characterize the slow dynamics in¬ 
duced by a perturbative single-photon loss on this protected subspace, by 
performing adiabatic elimination in the Lindbladian dynamics. 


1 Introduction 

The harmonic oscillator is a standard quantum system. It features coherent states, 
equivalent of classical harmonic oscillator amplitudes, whose coherent superposi¬ 
tions also called “cat states” feature genuinely quantum properties with no classical 
equivalent. In a recent paper [8] , our collaborators take advantage of this fact to 
propose an implementation of logical quantum bits as “cat states”, with potential 
to realize universal quantum computation using standard technological elements 
as quantum gates. Their key contribution, besides the insight that the cat states 
are inherently insensitive to part of the typical perturbations, is the design of an 
“engineered reservoir” that stabilizes a subspace of such states in open loop. 

The contribution of the present paper is to precisely establish, from the Lind¬ 
blad master equation (2), the stabilization properties of this scheme both, 

• in ideal situations with e = 0, by proving in theorem 1 global convergence 
of the infinite-dimensional nominal model towards the target “protected sub¬ 
space”, 
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• in presence of a small but dominant decoherence source with e > 0, by 
establishing the approximate slow dynamics as a reduced Lindblad master 
equation (15). 

The first goal builds on a typical Lyapunov-LaSalle strategy, with additional care 
needed due to the infinite dimension. For the second goal, we resort to a sepa¬ 
ration of the quantum dynamics into fast and slow components. We then apply 
an adiabatic elimination of the fast system to deduce a good approximation of the 
dynamics on a slow manifold, with the state remaining e-close to the protected sub¬ 
space. Studying such perturbations is standard for quantum Hamiltonian systems, 
where regular perturbation theory can be routinely applied [11], but the Lindbla- 
dian case with singular perturbations has attracted much less attention. In [9] and 
similarly [10] singular perturbations up to second order are applied to a system 
with N ground states and eliminating relaxing excited states. In [1, 13] specific 
atom optics dynamics and an ancilla-mediated feedback are investigated with the 
standard approach of [2]. In [6] the so-called Schrieffer-Wolff formalism is gen¬ 
eralized to Lindbladian dynamics; its basic form requires inversion of the nominal 
dynamics operator, which is not too practical and which we circumvent here for 
the derivation of the reduced slow master equation (15). 

The paper is organized as follows. Section II describes the mathematical model 
of the dynamics to be studied. Section III provides the global convergence proof for 
an idealized model, and Section IV analyzes precisely how this allows to counter 
the dominant external perturbation, which is single-photon loss. We pedagogically 
present the corresponding slow/fast perturbative argument (Section IV.B) as the 
translation to quantum notation of the standard dynamical systems approach, re¬ 
called in Section IV.A. Finally, simulations illustrate the validity of our analysis in 
Section V. 

2 Driven dissipative pairwise photon process 

The underlying space of the quantum harmonic oscillator is a Hilbert space “K of 
infinite dimension spanned by the Fock states {|n)|„6N- The annihilation operator a 
is defined by a\n) = ^ln\n - 1) for any n>\ and a|0) = 0. Its Hermitian conjugate 
a\ verifies a^\n) = ^n + \\n + I>, for any n > 0; . We denote N - a^a the photon 
number operator, satisfying A^|?i) = n\n), for any n > 0. For any a € C, a coherent 
state \a) € TF is characterized by a\a) - a\a} and 



The coherent states \a) are viewed as the quantum model for a “classical state” of 
complex amplitude a. Indeed, whereas the classical undamped harmonic oscillator 
j^x = -ojp, j^p = cox with X, p € R has solutions x{t) - the quantum 

harmonic oscillator ^|(A) = -ic){N + //2)|i/^) features solutions |i/^(0) = |Q'(0) - 




up to an irrelevant phase. As is customary, in the rest of the paper we 
describe the system in a frame rotating at the oscillator frequency co, for which the 
oscillator has stationary solutions \a). 

We call Schrddinger cat state or simply cat state the coherent superposition of 
two coherent states with opposite amplitudes, 




b)± I - a) 

y± 


( 1 ) 


where y- - is a normalization factor. If \a\ » 1 then we have 

w w V2. 

Throughout this paper, we denote by '7C^(‘K) the set of trace-class operators 
on “K, i.e. compact Hermitian operators on 'K whose eigenvalues {(Tk)keH satisfy 
Zfeo If/tl < +°°- This equipped with the trace-norm Tr(| • |) is a Banach 

space (see e.g. [12]). The set of density operators (quantum states) corresponds to 
elements of '7C^(‘K) that are non-negative and of trace 1. They are usually denoted 
by p. We denote also by the subspace of of operators whose range 

is included in a vector space spanned by a finite number of Fock states: 


fn,n'\n)(n\ : e CJnp = 

I finite 

We consider the quantum harmonic oscillator interacting with its environment 
as described in [8]. An external coherent driving field of amplifude u (assumed 
here real and sfricfly positive wifhouf loss of generalify) is applied such fhaf fhe 
oscillafor can only exchange photons in pairs. Furfhermore, fhe quanfum system is 
builf such fhaf similarly fhe main dissipative process is a pairwise phofon loss wifh 
rate /<• > 0. Nevertheless, due to physical constraints, even if a single photon loss 
process is much less frequent than the previous one, it can’t be totally omitted. For 
this reason, this term appears with a small coefficient 0 < e /c in the following 
Lindblad master equation which governs the system dynamics: 

^ - a^,p] + K2a^(p) + e2a(p) 

at 

where [•, •] stands for the commutator and where, for any linear operator A on 'K, 
the super-operator is given by 

£a(p) = ApA^ - (A^Ap -I- pA"^A)l2 . 

Noting a - 2uI k and L = one can reformulate the previous equation as : 

d 

—p = k2l{p) + e2a{p) 
dt 

It is shown in [8] that, for 6 = 0, the two-dimensional Hilbert space 

'Ha = span{|a>, | - a>} 


( 2 ) 




is a decoherence-free space: for e = 0, any density operator p with support included 
in "Ha is a steady state, i.e., fiitjd) = 0. 

We will not investigate here the well-posedness of (2) and the associated strongly 
continuous semigroup of linear contraction on TC'CK). Such issues can be inves¬ 
tigated via theorem 3.1 of [3] ensuring the existence of minimal solutions since 
for any a,K,e > 0, the operator -kL^L - ea^a is the infinitesimal generator of a 
strongly continuous one parameter contraction semigroup on TY. This means that, 
in the sequel, we will always assume that the Cauchy problem (2) with an initial 
condition po e positive semidefinite and of trace one, admits a solution in 

defined for any t > 0. Moreover, fhis minimal solufion will remain posifive 
semidefinife. We will invesfigafe here fhe time asympfofic regime for e = 0 and 
fhen for 0 < e /c. 


3 Convergence of (2) for 6 = 0 

Lemma 1. For any quantum state p in'K^Fl) we have Tr{^L2L{p)L"'^ < -2Tr{LpL^^. 

Proof. From Tr (L£t(p)L^^ = Tr^LpL^[L^,L]^and [L^,L] = [(a^)^,a^] = -4A^- 
21, we have 

Tr(L£L(p)L^) = -2Tr (LpL’f(2A^ -t /)) . 

We conclude since N and LpLf are posifive semidefinife. □ 

Modulo argumenfs proper fo fhe infinife-dimensional seffing, fhis basically 
means fhaf V{p) = Tr {hpLf^ is an exponenfial Lyapunov funcfion for fhe sysfem 
(2) wifh e = 0. 

Lemma 2. For any v > 1 there exists p > 0 such that, for any quantum state p in 
'KfFl) we have 

v+1 

rr(£i(p)iV0<-v(rr(piV0) +p 

Proof Denofe by fhe adjoinf of £/,: £^(A) = Lf AL - {if LA + ALf L)I2 for 
any Hermifian operator A on FI. Compufafions relying on fhe identify af{N) - 
f{N + I)a for any funcfion /, yield 

^liRN)) = -N{N-l){f{N)-f{N-2I)) 

+ ^a\f{N) - f{N - 21)) + ^{RN) - RN - 2ma^f . 

Since Tr(£i(p)/(iV)) = Tr (pfi^C/CiV))), we have 

= -Tr (pN{N - l)(/(iV) - RN - 21))) 

+ "^i:r(pa^{f{N) - RN -21))] 

+ ^^Tv(p{RN)-RN-2ma^f). 


Tr(£i(p)/(iV)) 


Take f{x) = and define g{x) - f{x) - f{x - 2) for x > 2; g{x) - f{x) for 
2 > X > 0; g{x) = /(O) for x < 0. From a^g(N) = a yjg{N - I)a yjg{N) and 
Cauchy Schwartz inequality 


|Tr (pa ViV)) I - 


< 


Tr |( ^|pa yjgiN - /))(a x/g{N) j 
^Tr (pg{N - 2I){N + /)) Tr (pg{N)N). 


Since g{x - 2)(x + 1) < g(x){x + 3) for x > 0, we have 


|Tr (pa2g(iV)) I < Tr (pgiNXN + 31)). 


Thus 

iTr (p(a2g(iV) + g(iV)(a1^)2)) < Tr ipg{N){N + 31)) 
and 

Tr (£l(p)A^ 0 ^ Tr (pg(iV) (-iV^ + (^2 + + 3 ^ 2 ^ j 

Since (x^ - (x - 2)^) (-x^ + {a^ + l)x + 3a^^ is equivalent to -2vx''''‘' for large x, 
there exists p > 0 such that for all x > 0, 

((x + ly - x^) (-x^ + (a^ + l)x + 3a^^ < + p. 

Finally we get 

Tr{2L{p)N'') < -yTr(pAt''+') +p 

V 

Since x^x is concave, ^Tr > Tr(pA^’'). □ 


Theorem 1, Consider a trajectory [ 0 ,+oo [9 t 1 -^ p{t) € 'K^yH) of the master 
equation (2) with e = 0, k> 0 and a > 0. The following statements hold true: 


1. Take v > 1. Then there exists y > 0 such that, for any initial quantum state 
p(0) = po satisfying Tr {poN'') < + 00 , we have for all t > Q, Tr {p{t)N^) < 
max {y,Tr(poN'')y 


2. Assume that po € ’K^ifH). Then, there exists a quantum state p with support 


in Tia such that, for any y > 0, lim,H^+oo 7>( (p{t) - p)N 2 


- 0 . 


The limit p depends on po- It can be derived from po with the four Hermi- 
tian, bounded and independent operators that are in the kernel of the adjoint super¬ 
operator and given in [ 8 ]. 











Proof. The first statement is a direct consequence of Lemma 2 and of the fact that 
quantum states element of are a dense subset of the quantum states cr of 

(TT) with Tr {aN'') finite. Tr remains bounded since 

-f-Tr(piV'') - kTt(2l(p)N'') < -vkTt + ku. 
at 


Thus for Tr(pAf'') >A = iTr(pN'') < 0. 


The second statement exploits the first one. We can assume v >2. Take v' > v. 


Denote by TClCK) the supspace of trace-class operators cr such that Tr 


is finite. The space 'KfifH) with the norm \\(t\\v> = Tr(|cr|) -i- Tr( 


iV2 criV2 


iV2 crN 2 


IS a 


Banach space. From the first statement, we know that po being an element of 
'Kl.ifH), p{f) remains always in 'K], {!-(). Since v' > v, the injection of 
into TCjCK) is compact: {p{t) | t > 0) is precompact in Denote by p € 

TCjCK) an adherent point of p(t) for t tending towards infinity. Since y > 2, p 
and p belong to the domain of £/,. Lemma 1 implies Tr^LpL^^ = 0, i.e., the 
support of p is contained in the kernel of L, which coincides with 9^a■ Moreover, 
the semigroup associated to the Lindblad master equation is a contraction for the 
trace distance: for two trajectories pi(0 and p 2 (^)^ t ^ Tr(|pi(t) -p2(0l) is a non¬ 
increasing function. Thus t i-^ Tr (|p(f) - p|) is non-increasing since p is a steady 
state. Consequently the adherent point p is unique: p(t) converges towards p in 

-x-.H-K). □ 


4 Reduced slow dynamics of (2) 

We have proved in the previous section that the system converges toward the deco¬ 
herence free subspace 9^a when we neglect the photon loss channel (e = 0). When 
0 < e <s 1, the center manifold theorem allows us to separate the system into fast 
and slow dynamics. The fast dynamics makes the system globally converge to a 
subspace close to ‘Kq.. The slow dynamics approximate the behavior of “protected 
states” |cj), |c“). The present section is aimed at characterizing these dynamics to 
the first order in e. 

The fast/slow dynamics reduction from nonlinear systems theory, also known 
as singular perturbation theory, is useful despite the linearity of Lindbladian dy¬ 
namics, because the very high dimension and often high degeneracy of open quan¬ 
tum systems makes matrix diagonalization impractical to apply. 

For the sake of clarity, we first particularize the fast/slow dynamics reduction 
theory to linear systems of finite dimension with the standard notations. We pro¬ 
pose then a reduction procedure via an adapted duality viewpoint with the charac¬ 
terization based on (9) and (10) here below. We apply then this characterization 
to the quantum system (2) and show that it facilitates the computations up to first 
order with respect to e.g. [6]. 






4.1 Reducing a linear system to its slow dynamics 

We here review the theory of geometric singular perturbation, which was mainly 
developed by Fenichel in [4] and surveyed by Jones in [5], in a linear context. 
Consider the nominal linear system j^x = Ax with x = {x\,X 2 ) e R™ x R" and 
converging globally to the m-dimensional subspace .S = {x £ R™’*’” : X2 £ R” = 0). 
To this nominal dynamics we add an arbitrary perturbation matrix B of order e <sc 1. 
In matrix notation, the dynamics can be written in block form which yields: 

j^Xi = AiX 2 +e{ByX2 +BqXi) (3) 

J^X 2 = A 2 X 2 + e{B 2 X 2 + B 2 ,Xi) . 

The assumption that S is globally exponentially stable for e = 0 corresponds to A 2 
having all eigenvalues with strictly negative real parts, thus it is invertible. Hence 
we can define the regular change of variables 

Xi - Xi - A\A2^X2 
X2 = X2 

which yields dynamics in Tikhonov normal form 

= e{{BQ - A\A2^B2 ,)x\ 

+{Bq + B\- AiA2^{B2 + B2,))X2^ 

= ef{xi,X2). 

j^X2 = A 2 X 2 + e{B2X\ + {B 2 + B3AiA2^)x2) 

= g{xi,X 2 ,e). 

The Tikhonov conditions for reducing the system by singular perturbations is that 
the first (slow) subsystem has eigenvalues going down as e, while the second (fast) 
subsystem has eigenvalues bounded away from zero for e = 0. These conditions 
are satisfied above. The Tikhonov fheorem fhen allows fhe following reducfion. 

Proposition 1. The trajectories of the full system (5), (6) (with initial conditions 
satisfying g = 0) remain e-close over at least a time of order 1/e, to the trajecto¬ 
ries of the system restricted to the “slow submanifold” g{x\,X 2 , 0 ) = 0 and where 
dynamics are given by replacing X 2 in f{xi,X 2 ) by the solution of g(x\,X 2 , 0) = 0. 

In our linear case, fhe slow manifold comes down fo X2 = 0 and fhe slow 
dynamics frivially reduce fo fhe firsl term in (5). Transforming back fo fhe original 
coordinafes fhe slow manifold corresponds jusf fo X2 = 0 and fhe dynamics are 

-e(Bo-A 1A2^63)^1- (7) 

The second term rellecfs fhe influence of fhe fasl X 2 dynamics on fhe slow variable 
x \: by blindly setting X2 = 0 in fhe original system (3) and neglecting ifs second 


(4) 

(5) 

( 6 ) 


line, we would miss this term and get an incorrect approximation. 


Computing the corrective term by explicit inversion of A2 can be 

a tedious task when the fast subsystem has a large dimension (in our quantum 
case, X 2 would rigorously be of infinite dimension). However if first integrals of 
the system with e = 0 are known, these can facilitate the computations via the 
following dual viewpoint. 

Consider a linear functional = (p^, p^) e which is conserved by 

X = Ax, i.e. satisfying p - A^p = 0 or equivalently p^AiX 2 + P 2 A 2 X 2 = 0 for all 
X 2 € R”. (The notation denotes matrix transpose.) Again using invertibility of 
A2, we see that p^ actually satisfies 

p[AiA2 V + ply ^OVy eR" . (8) 

Knowing m linearly independenf funclionals {p^{k))k=i,...,m salisfying (8) is sufh- 
cienf fo fully characterize the corrective term in (7): we have 

= 6(Bo + Q)xi (9) 

with Q defined by fhe sef of linear equations: 

pl{k)Q = pl{k)B 2 , k = \,2,...,m. (10) 


4.2 Quantum system (2) with 0 < e /r 

We know apply fhe same procedure fo our quanfum system. The nominal x = Ax 
corresponds to p = ^hip)- We have shown in Section 3 that: 

• S corresponds to a four-dimensional real subspace of Hermitian operators 
spanned by |c+>(c+|, |c“)(c“|, |c+><c“| -1- |c“>(c+|, and /(|c+>(c“| - |c“)<c+|). 
These correspond to the m = 4 coordinates of xi. 

• The subspace S is globally asymptotically stable under p = Ql{p)- This 
corresponds to the invertibility condition on A2 independently of e. 

We therefore introduce the projector 

Pc = \ct}{4\ + K)(c~a\ (11) 

such that 

p, = PcpPc (12) 

corresponds to the “slow” xi space of the previous section. The projection onto 
“X 2 space” is given by 

Pf =p- PcpPc ■ 

Our goal is to compute the evolution of p^, which is the equivalent of (7). For 
this we take advantage of the dual formulation (9), (10). The following procedure 


can in principle be applied to any perturbative dynamics, we here focus on as a 
physically relevant case. 

The perturbative dynamics on the slow manifold features a first component, 
corresponding to Bq, obtained simply by projection onto the slow manifold. The 
identities a\c-) = quickly yield its explicit expression: 

PA(Ps)Pc - a^^xiPs) (13) 

where X - — Ic+Xc"! + — Ic-Xc+I. (14) 

7- 7+ 

To compute the corrective term by duality, using the equivalent of (9),(10), we 
need to identify m - 4 conserved functionals of the system, which are the fixed 
points of the dual nominal dynamics Fortunately, those invariants are 

known for the particular operator L: 

• One easily checks that = I the identity operator is in the kernel of any 

• The parity operator = (-1)“ “ is in the kernel of because photons are 

exchanged by pairs. 

• The appendix of [8] gives two more operators and in terms of Bessel 
functions; one checks that they are linearly independent for finite a. 

We will also use the following key property of the conserved quantities, which is 
specific fo the structure of quantum Lindblad dynamics. 

Lemma 3. Any Hermitian operator ^ in ker(fip commutes with Pc the orthogonal 
projector onto j-(a- 

Proof From 2Lf^L - Lf Lf + and LPc = 0 = PcL^ we have L^L^Pc = 
0 = Pc^L^L. Thus the Hermitian operator A = Pc^ + ^Pc satisfies LALf - 
(L^LA + AL^L)I2, i.e. belongs to the kernel of £/,. The support of A is thus 
included in jia and thus [Pc, A] - 0. This implies that Pc^ - PcfPc = ^Pc- □ 


Considering projections on respective subspaces, the equivalent of the term 
of the linear perturbation is given by QaiPs) - Pc^a(Ps)Pc- The equivalent of equa¬ 
tion (10) characterizes the Hermitian operator Q with support on jPa as follows: 


TriPcePcQ) 


Tr({e - PcfPc) (fi.(p.) - PA(p.)Pc)) 

Tr[e(^a(p.)-PA(Ps)Pc)) 

iTr (r {(Pc - + p.a'^aiPc - /))) 

iTr 1^'' (a'^ap,{Pc - /) + (Pc - I)p,a^afj - 0 


for any v = a, b, c, d. From the first to the second line, Pc^Pc is readily dropped 
since Pc is a projector. For the next one we have to write out and see that 


PcUpsO^Pc - apsO.^ for the particular perturbation operator a. The last line follows 
by using Lemma 3 and PcPs - Pi- As a conclusion, we get that the corrective term 
is Q = 0 for our particular case. We can summarize these computations as follows. 

The trajectories of system (2) with initial conditions having supports in Tia 
remain e-close over at least a time of order 1/e, to the trajectories of the “slow 
variable” ps which is a linear combination o/|c|^>(c^|, |c“>(c“|, |cj)(c“| + |c“>(cj|, 
and /(|cj)(c“| - |c“)(c^|). The slow state ps follows the Lindblad dynamics: 

j^p, = ea^&xiPs) (15) 

where X - — |c'*')(c“| + — |c“)(c'''|. (16) 

r- r+ 


In applications considering |c^) and |c“) as canonical states |0), |1) of a logical 
qubit [8], the operator X corresponds to a bit-flip in the limit ^ > 1 of large 

coherent amplitude a and to a decoherence to the vacuum |0) in the limit of Fock 
states \cjf) = \n = 0), |c“) = \n = 1) when a = 0. For all other cases, the qubit 
dynamics (15),(16) corresponds on the canonical Bloch sphere to: 


d 

2 (yl - r-)" 



= 

-a - — X 

2ylyl 




2 (.yl^ylf 



d,y = 

2ylyl " 



A. - 

dt^ 

2 rt + y- i 
® 2 2^ 
y+y - '' 

rt-yf) 
rt+yf / 


This converges to x = 0 (slowly 

for ^ l),y ^ 0 , 

and z = 

= ^4 (which is - 0 for 

yj+yl 


- 1 ). 

r- ' 

5 Numerical simulations 

To illustrate the interest of this model reduction based on adiabatic elimination of 
the rapidly converging variables, we compare via numerical simulations trajecto¬ 
ries of the complete system and of the reduced one. To do so, we use a numerical 
scheme which preserves the positiveness for the Lindblad equation and similar to 
one used in [7]. We choose the following values of the parameters : a time-step of 
10~^, the decoherence strength k = 1, u = 1/2 and e = 0.01. With a = 1, the pop¬ 
ulation of photons for n > n^ax = 40 in the coherent state |a) is almost zero since 
less that Consequently, we truncate the infinite-dimensional Hilbert space 

to a numerical system space spanned by {|1), |2),.. . |40)) in the Fock basis. We 
denote by p the density matrix of the resulting system. It is worth stressing that the 
complete system is represented by an n^ax x n^ax matrix while the reduced system 







is represented by a 2 x 2 matrix (on the basis |cj>, |c“)). Thus the computation is 
much faster on the second one. 

We first take as initial condition the vacuum state, po = |0)(0|- The state of the 
reduced system, ps is then initialized at |c|^>(c^| because both po and |cj>(c^| are 
+ 1 eigenstates of the parity operator = (-1)" which is a conserved quantity 
(see Section 4.2 and [8]). To compare the trajectories of (2) initialized at |0)(0| and 
of (15) initialized at |c^)(c^|, we show in figure 1 fhe expecfafion values Tr{pcr^ 
and Tv{psO-^ of fhe operator ctj, - |cj)<c^| - |c“>(c“|, commonly denoted (cr^). 
Afler a fransifional regime of fypical durafion l/zr, one can see a sfrong similarify 
befween Tr (pcr^) and Tr (p^cr^) up fo a consfanf offsel. The value of fhis offsel is of 

order e. Furfhermore, we plof fhe fidelily F(p,ps) - t'" | befween p^ 

and p. For beffer readabilify, figure 2 shows fhe logarifhm of 1 minus fhe fidelify, 
i.e. of ifs deviafion from fhe ideal value 1. This deviafion quickly converges fo an 
order 10““^, corresponding to as expected. If fhen furlher decreases, incidenfally, 
as bofh sysfems converge towards fhe unique equilibrium of fhe slow dynamics. 

To emphasize fhe influence of y''' and y~ in (16), we add a simulafion wifh fhe 
same parameters buf wifh fhe following and same initial condition for fhe complefe 
and reduced system: 


p^ = \ (l4) + lca>) (<4l + <c«l) 

Figure 3 shows fhaf fhe expecfafion value of (Tx = + |c“)(c^| slowly de¬ 

creases over time, as expecfed from bif-flip dynamics. The slope of fhis decrease is 
approximated fo ~ 4% accuracy by fhe reduced dynamics. Moreover, figure 4 es- 
fablishes fhaf {cr^) does nof remain zero. This is due fo fhe facf fhaf, wifh y'*' > y~, 
equation (15) “promotes” fhe population of |c^)(c^| over fhe population of |c“)(c“|, 
unlike a pure bif-flip. 

The simulafions fhus confirm fhe validify of our approximafion of fhe complefe 
model by fhe reduced one. 

6 CONCLUSIONS 

We have rigorously proved convergence of a harmonic oscillafor Lindblad dynam¬ 
ics wifh fwo-phofon exchanges, fo a protected subspace. We have also esfablished 
fhe approximafe slow dynamics on fhis profecfed subspace when a fypical perfur- 
bafion is added, and illusfraled ifs validify in simulafions. The mefhods used for 
fhis particular example are applicable to general Lindbladian dynamics. 

The reducfion by singular perfurbafions and adiabafic eliminafion is of course 
applicable in general fo evaluate fhe remaining slow dynamics in quanfum sysfems 
wifh (engineered) profecfed subspaces. Extension fo k-phofon processes wifh 
k > 2 can be addressed in fhe same way. The facf fhaf fhe slow variable sfill fol¬ 
lows a Lindbladian master equafion may nof be surprising buf remains fo be proved 




Figure 1: Comparison of Tr (cr^p) and Tr (cr^s) for p solution of the complete 
system (2) with e - a = 1 and vacuum initial condition (truncation up to 
40 photons), and for ps solution of the reduced system (15) with initial condition 

14 X 41 . 



Figure 2: logjQ(l - F) where F is the fidelity between ps and p for the simulations 
of figure 1 . 














Figure 3: Comparison of Tx{(TxP) and Tr{o-xPs) for p solution of the complete 
system (2) with e - a = 1 (truncation up to 40 photons), and for ps solu¬ 
tion of the reduced system (15), with the same initial condition p(0) = Pi(0) = 
5 (|c+>-i-|c“>)((c+|-i-(c“|). 



Figure 4: Comparison of Tr (cr^p) and Tr {(r^ps) where p and ps correspond to 
simulations of figure 3. 















in the general case. The fact that the dynamics reduces to the orthogonal projec¬ 
tion of the Lindbladian onto the protected subspace (i.e. Bq without any correction 
due to B 3 , in the terms of Section 4.1) for the case examined here would be in 
agreement with the physicists’ “quantum Zeno” viewpoint. However, under which 
formulation this viewpoint should be applied in the general case also remains to be 
rigorously characterized. 
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